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AM  APPROXIMATE  THEORY  FOR 
TURBULENT  WAKES  BEHIND  HYPERSONIC  BLUNT  BODIES 


W.  W.  Short 


ABSTRACT 

The  growth  of  the  turbulent  core  in  the  constant  pressure  wake, 
of  a  blunt  hypersonic  body  is  predicted  by  means  of  a  mixing  length 
theory.  Near  the  body  where  the  core  is  smaller  than  the  enthalpy 
wake  the  constant  density  theory  of  Townsend  is  applied.  A  variable 
density  solution  is  obtained  for  the  region  in  which  the  core  and  en¬ 
thalpy  wake  are  of  the  same  size.  The  two  solutions  are  fitted  to 
predict  velocities,  temperatures,  and  concentrations  in  the  constant 
pressure  region  for  several  thousand  body  diameters  behind  the 
missile. 

I.  INTRODUCTION 

The  environment  around  hypersonic  bodies  in  the  atmosphere  is  not  well 
understood  especially  in  the  base  and  wake  regions.  In  spite  of  the  uncertainty 
associated  with  the  chemical  processes  occurring  in  high  temperature  air,  theo¬ 
retical  analyses  of  the  shock  layer  and  boundary  layer  have  met  with  a  great  deal 
of  success;  whereas,  analyses  of  the  base  and  wake  flow  have  not  been  particularly 
rewarding.  This  is  unfortunate  when  attempting  to  detect  and  discriminate  between 
bodies  entering  the  earth's  atmosphere  because  the  wake  comprises  the  largest 
portion  of  the  observable  flow  field.  In  the  present  analysis,  a  simple  model  of 
the  turbulent  wake  has  been  used  to  calculate  the  properties  in  the  wake.  The 
effect  of  a  few  foreign  materials  in  the  wake  is  discussed  because  they  may  have 
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a  profound  effect  on  the  electromagnetic,  radiative,  and  transport  properties  of 
the  flow.  A  foreign  material  (or  additives  as  they  are  referred  to)  is  any  chemi¬ 
cal  substance  not  normally  found  in  the  atmosphere.  An  understanding  of  the 
effect  of  additives  on  the  environment  around  hypersonic  vehicles  may  be  neces¬ 
sary  to  solve  the  problems  of  detecting,  discriminating,  and  destroying  re-entering 
ballistic  missiles. 

1  2 

Feldman  and  Goulard  have  obtained  solutions  for  the  physical  properties 
in  a  laminar  trail  behind  a  missile.  Laminar  wakes  are  very  long;  the  length 
being  proportional  to  the  frontal  area  of  the  missile.  Feldman's  calculations  show 
that  at  an  altitude  of  60, 000  ft.  and  a  velocity  of  17, 500  ft/sec  the  temperature 
1,000  km  behind  a  1-m-diam  missile  Is  1500 **0,  or  about  one-half  the  wake  tem¬ 
perature  near  the  missile.  Laminar  solutions  also  indicate  that  the  width  of  a 
hypersonic  wake  is  nearly  constant.  Townsend  developed  a  mixing  length  theory 
which  predicts  velocities  in  the  turbulent  wake  behind  a  body  passing  through  an 
incompressible  fluid.  In  this  case  the  length  of  the  trail  is  proportional  to  the 
diameter  of  the  missile,  not  its  area,  and  the  width  of  the  turbulent  zone  increases 

as  the  cube  root  of  distance  behind  tlie  body.  Tovmsend's  incompressible  theory 
4  0 

was  modified  ’  and  applied  to  hypersonic  wakes.  The  present  work  is  based  on 
this  modification. 

The  problems  of  transition  from  laminar  to  turbulent  flow  in  the  flow  field 
are  not  considered.  It  is  assumed  the  wake  is  turbulent  over  its  entire  length. 

The  wake  is  not,  however,  turbulent  over  its  entire  width  at  all  points. 


n.  FLOW  FIELD 

Before  presenting  the  theory  of  turbulent  diffusion  in  a  hypersonic  wake  the 
properties  of  the  wake  will  be  discussed  qualitatively.  The  general  properties  of 
a  hypersonic  wake  are  outlined  in  Figure  1.  The  flow  behind  a  missile  can  be 
divided  in  the  axial  direction  into  three  regions;  the  base  region,  the  expansion 
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FIGURE  1  TURBULENT  WAKE  BEHIND  A  HYPERSONIC  BODY 


wake,  and  Uie  constant  pressure  wake.  The  constant  pressure  region  is  often 
termed  the  far  wake.  The  wake  can  bo  further  divided  in  the  radial  direction  into 
the  viscous  and  inviscid  flow.  The  turbulent  core  constitutes  the  region  in  which 
viscous  forces  are  important.  Surrounding  the  core  is  the  inviscid  region  which 
is  caused  mainly  by  Uie  lx)w  shock.  It  is  the  groulh  of  the  turbulent  zone  into  the 
inviscid  flow  and  the  consequent  diffusion  of  the  high  enthalpy  trail  which  will  be 
discussed  here. 

The  Inviscid  flow  behind  a  blunt  lx)dy  can  be  characterized  by  a  maximum 

velocity,  U  ,  which  is  about  one-fifth  of  the  missile  velocity,  a  maximum  tem- 
0 

perature,  T  ,  which  is  about  one-half  of  the  stagnation  temperature,  and  a  width, 
0 

L  ,  which  is  several  times  that  of  the  body.  Figure  2  shows  a  cross  section  of 
0 

the  viscous  and  inviscid  portions  of  the  wake  behind  a  spherical  body  having  a 
velocltj'  of  17,500  ft/sec.  The  large  variations  in  velocity,  temperature,  and 
concentration  near  the  a.vls  originate  in  the  boundary  layer  at  the  surface  of  the 
missile.  The  boundary  layer  on  the  body  will  flow  Into  the  turbulent  core  at  the 
axis  of  the  wake.  Temperatures  in  the  turbulent  core  are  sensitive  to  the  type  of 
surface  material  on  the  body.  If  the  body  surface  were  adiabatic  the  gas  tempera¬ 
ture  In  the  base  region  and  at  the  a-\ls  of  the  core  immediately  downstream  of  the 
base  would  be  verj'  close  to  the  stagnation  temperature.  In  practical  cases,  the 
surface  material  absorbs  energy  from  the  boundary  layer,  and  the  gases  radiate 
energy  to  the  surroundings;  therefore,  the  temperature  at  the  axis  of  the  turbulent 
core  may  be  closer  to  that  of  the  missile  surface  than  the  stagnation  temperature. 
All  material  which  flows  into  the  boundary  layer  on  the  missile  will  be  deposited 
in  the  turbulent  core.  Therefore,  material  ablated  from  the  missile  surface  will 
be  deposited  in  the  core  as  shown  by  the  concentration  curve  in  Figure  2.  Thus, 
only  the  properties  of  the  boundary  layer  and  core  of  the  wake  will  be  affected  by 
ablation.  However,  some  properties,  the  electron  concentration  for  example, 
may  be  very  sensitive  to  small  amounts  of  foreign  materials  present.  The  dimen¬ 
sions  and  thermodynamic  properties  of  the  core  are,  therefore,  of  major  interest. 
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Dimensionless  Radius,  r/R 

FIGURE  2  VELOCITY,  TEMPERATURE,  AND  CONCENTRATION 
BEHIND  A  HYPERSONIC  BLUNT  BODY 
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ra.  ASSUMPTIONS 


It  is  assumed  in  this  analysis  that  the  core  of  the  wake  is  completely  turbu¬ 
lent.  This  assumption  is  based  on  pictures  obtained  In  hypersonic  ranges  at  1  atm. 
The  possibility  of  transitions  from  laminar  to  turbulent  flow  occurring  In  the  wake 
behind  a  missile  at  high  altitudes  is  neglected.  Pressure  variations  in  the  expansion 
wake  are  not  included  in  tlie  luialysis.  The  results  of  radiation  measurements  made 
by  Georgiev  in  a  hypersonic  range  indicate  that  the  point  at  which  the  pressure  is 
nearly  equal  to  ambient  pressure  is  appro.xlmately  40  diam  behind  the  missile. 
Errors  contributed  by  the  a.ssumptlon  of  constant  pressure  in  the  first  40  diam 
will  not  greatly  affect  the  dimensions  of  the  total  wake  which  is  several  thousand 
diameters  in  length.  The  Increase  in  temperature  by  viscous  dissipation  in  wakes 
was  found  by  Feldman^  to  be  small  and  is  neglected.  Temperatures  in  a  wake  are 

in  general  less  than  3000 °K  for  re-entry  speeds;  therefore,  the  wake  was  treated 
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as  a  perfect  gas.  The  thermod)mamlc  properties  of  high  temperature  air,  however, 
were  used  to  obtain  the  initial  properties  of  the  wake. 

The  exact  nature  of  the  turbulent  core  will  depend  upon  the  geometry  and  the 
flight  conditions  of  the  body.  Because  of  the  wide  variety  of  temperature  profiles 
possible  in  the  turbulent  core,  it  has  been  assumed  to  be  Initially  isothermal  in 
this  calculation.  Therefore,  the  turbulent  core  analysed  here  will  consist  of  a 
sharp  velocity  profile  as  depicted  in  Figure  2  in  a  uniform  temperature  and  pres¬ 
sure  field.  The  early  growth  of  the  turbulent  core,  in  this  case,  can  be  estimated 

3 

by  the  incompressible  theory  developed  by  Townsend.  It  is  assumed  that  this 
analysis  is  valid  to  the  point  at  which  the  width  of  the  turbulent  core  is  equal  to 
that  of  the  inviscid  wake;  beyond  this  point  a  different  analytical  solution  will  be 
obtained. 

The  assumptions  used  are  summarized  by  the  following  list: 

1.  The  fluid  is  a  perfect  gas; 
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2.  The  wake  is  axl-symmetric; 

3.  The  core  is  completely  turbulent; 

4.  Constant  pressure, 

5.  Negligible  radial  velocit>'  components; 

G.  Steady  sta;e; 

7.  Negligible  viscous  dissipation; 

8.  Prandtl  mixing  length  theory'; 

C  € 

9.  Eddy  Prandtl  number  Is  unity,  ^  =  1 


IV.  ANALYSIS 


Townsend  has  showm  by  means  of  a  mixing  length  theory  that  when  self- 
preser\‘ing  flow  is  established  the  velocity  distribution  is  given  by 
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u  =  u  exp 
a 


1  r_ 

2  2  ! 


(1) 


The  nomenclature  is  listed  on  p.  17 .  The  lower  case  u  represents  the  velocities 

in  the  core  relative  to  the  maximum  velocity  of  the  inviscid  flow  at  the  same 

point,  U,  assuming  the  core  did  not  exist.  The  total  velocity  at  any  point  Is, 

therefore,  u  +  U.  The  width  of  the  core,  I ,  and  the  velocity  at  the  axis,  u  ,  are 

a 

given  by 
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Equations  (2)  and  (3)  Include  the  drag,  F,  which  in  the  hypersonic  case  is  the  drag 
experienced  by  the  core  which  is  the  skin  friction  and  base  drag  but  not  the  wave 


7 


drag.  Therefore,  F  can  lie  e.vpressed  in  terms  of  the  IniUal  width  and  velocity  of 
the  core  by  using  the  following  momentum  integral 


F  = 


(4) 


A  solution  to  this  integral  equation  can  be  obtained  by  assuming  a  Gaussian  velo¬ 
city  dlstribuUon  as  used  by  Towmsend  in  Equation  (1).  The  drag  is 

F  -  r  (5) 


The  core  width  iuid  maximum  velocitj'  can  now  be  e.vpressed  by 


.3 

i 


3 ;  “  (X  -  X  ) 
o  0 


2  R 

t 


and 


„2  3  2 
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u  =  - 

a  2 
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(V) 


Equations  (1),  (6),  and  (7)  provide  an  approximate  method  for  predicting  velocities 
in  a  turbulent  wake  to  the  point  where  the  width  of  the  turbulent  core  has  grown  to 
equal  that  of  the  inviscid  wake.  The  turbulent  core  has  nearly  disappeared  from 
the  velocity  profile  because  the  momentum  of  the  core  is  a  small  fraction  of  the 
total  wake  momentum. 


A  solution  will  now  be  obtained  for  the  decay  of  the  fully  turbulent  compres¬ 
sible  wake.  This  solution  and  that  given  above  can  be  joined  to  give  a  complete 
description  of  the  wake.  By  using  the  assumption  listed  earlier  the  momentum 
equation  can  be  simplified  to  give 
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and  the  energy  equation  Is 


Sr 


U  ^ 
dx 


+  V 


1  d 


(9) 


The  distances  x  and  r  can  be  considered  dimensionless  since  each  term  contains 
X  ^  or  r  \  When  using  the  concept  of  eddy  viscosity  and  eddy  conductivity  the 
shear  and  energj'  flux  are  defined  by 


and 
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9  r 
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(10) 


(11) 


Following  the  mixing  length  theory  used  by  Townsend  it  is  assumed  that  the  eddy 
viscosity  is  proportional  to  the  wake  width  and  the  greatest  velocity  difference 
across  the  wake. 


The  following  relation  will  be  used 


e 

P 


(12) 


Where  R  is  the  same  unknown  constant  obtained  for  the  incompressible  case. 
Using  this  relation,  the  fact  that  h  =  C^T  for  a  perfect  gas,  and  assumption  5 
Equations  (8)  and  (9)  can  be  simplified  to  give 

9U  9  /  9U\ 

<»9x  Tp  dT\  ^  dT  J  (13) 


and 


U 


9T 

9x 


(14) 


The  velocity  term  on  the  left-hand  sides  of  the  momentum  and  energy  equations 

should  be  U;  however,  it  was  assumed  to  be  U  because  the  maximum  velocity  at 

00 

the  axis  of  the  fully  turbulent  wake  is  within  20%  of  the  free  stream  velocity. 
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The  solution  of  Equations  (13)  and  (14)  will  be  similar  provided  the  eddy 
PrandU  number  is  unity  and  the  Initial  conditions  are  similar.  If  these  two  con¬ 
ditions  are  true  the  temperatures  and  velocities  In  the  wake  can  be  related  In  the 
following  manner: 


U  =  a  (T  -  TJ  ^15) 

Substitution  of  Equation  (15)  into  Equation  (12)  gives 

a  L  (T  -  T  ) 
a  of 

V  "  (16) 

Comblnliig  Equatloa  (16)  and  the  perfect  gas  law  with  Equation  (14)  produces  the 
following  differential  equation  In  terms  of  temperature 

U 

t  °°  9T  _  T  3T\ 

~  VT  ary 

This  equation  can  be  solved  by  separation  of  variables  for  the  case  when  T  =  0. 

00 

This  case  Is  of  interest  because  the  ambient  temperature  is  very  low  (T  ~  200 “K) 
compared  with  temperatures  In  the  wake  of  re-entering  ballistic  missiles 
(T^  -  3000“K).  Assuming  that  T  =  f^  (x)  f^  (r)  gives  the  folloNvlng  relations: 


and 


(18) 


f2  =  be 


2aJ 


(19) 


It  is  believed  that  the  mixing  length  is  directly  related  to  the  width  of  the  wake. 
By  defining  the  mixing  length  in  compressible  flow,  L,  as 


00 
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the  variable  A(x)  can  be  evaluated  to  be 


A(x)  =  -  (X  +  c) 
Qt 


(21) 


where  c  is  a  constant.  Substitution  of  this  relation  into  Equation  (18)  gives 
an  equation  for  f^  in  terms  of  x.  The  temperature  is  approximated  by 


a  R  U 
t  “ 


(X  +  c) 


(22) 


when  «  T^.  Equation  (22)  only  satisfies  the  original  differential  equation 
(Equation  (17)  J  exactly  when  =  0.  According  to  Equation  (16)  when 

^erefore,  the  concept  of  eddy  properties  becomes  invalid. 

The  eddy  viscosity  can  be  greater  than  the  molecular  viscosity,  r  ,  but  it  can¬ 
not  be  less.  When  v  the  momentum  and  energy  equations  for  laminar 
flow  must  be  employed.  At  distances  greater  than  1000  diameters  the  constant 
width  wake  predicted  by  Equation  (22)  has  cooled  and  begins  to  grow  again; 
therefore,  this  analysis  cannot  be  extended  beyond  this  point. 

Wake  temperature  obtained  by  means  of  Equations  (1)  and  (22)  will 
always  follow  a  Gaussian  distribution.  This  would  require  that  the  initial 
co.idition  have  the  same  form.  The  radial  temperature  distribution  at  x  =  0, 
the  point  where  the  pressure  equals  the  ambient  pressure,  has  been  computed 
by  assuming  an  elliptical  shaped  shock  wave  around  the  nose  of  the  missile 

g 

as  suggested  by  Van  Dyke.  If  it  is  assumed  that  the  only  irreversible  pro¬ 
cesses  occurs  at  the  shock  wave,  the  thermodynamic  state  of  the  gas  at  x  =  0 
can  be  computed  by  following  an  isentropic  expansion  to  atmospheric  pressure 
along  each  streamline.  The  velocity  and  position  of  each  streamline  which 
passes  through  the  shock  wave  can  be  located  at  x  =  0  by  employing  the  Bernoulli 
and  continuity  equations.  The  result  of  such  a  calculation  is  shown  as  curve 
X  =  0  in  Figure  3,  which  gives  the  properties  of  the  inviscid  flow.  The 
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TEMPERATURES  IN  A  TURBULENT  WAKE 


properties  of  the  turbulent  core  were  obtained  by  assuming  u  =  U  -  u  and 

Its  momentum  equal  that  of  the  boundary  layer  on  the  body  aUts  p^lnt  of 
separation. 

The  constants  a  and  c  in  Equation  (22)  were  evaluated  by  fitting  the 
initial  temperature  profile.  It  is  obvious  in  Figure  3  that  the  curve  labeled 
X  =  0  does  not  follow  a  Gaussian.  A  finite  difference  solution  to  Equation  (17) 
using  the  conditions  at  x  -  0  was  carried  out  using  a  digital  computer  to  deter 
mine  how  quickly  this  curve  approached  a  Gaussian  distribution.  The  results 
of  this  numerical  solution  showed  the  initial  conditions  decay  to  a  curve  which 
can  be  fitted  by  Equation  (22)  within  20  diam.  This  decay  does  not  begin  until 
the  core  has  grown  to  the  width  of  the  inviscid  wake.  L.  A  constant.  X  . 
which  18  the  length  of  the  wake  required  for  the  core  to  grow  to  width  L^is 
added  to  the  dimensions  shown  in  Figure  3. 

Townsend  computed  theoretically  that  the  constant  has  a  value  between 
14  and  21.  A  value  of  14  was  used  in  the  present  calculations.  Results  ob¬ 
tained  by  means  of  Equations  (1).  (6).  (7).  and  (22)  are  shown  in  Figures  4 

and  5  for  a  missile  velocity  of  17.500  ft/sec.  an  altitude  of  60.  000  ft.  and 

an  initial  core  diameter  of  1/5  the  missile  diameter. 


V.  CONCLUSIONS 

The  present  analysis  is  approximate  and  subject  to  all  the  short  comings 
of  a  mbting  length  theory.  Several  sources  of  uncertainty  are  outlined  here. 

The  turbulent  dlftuslvlty  Is  assumed  to  be  a  function  of  the  wake  width 
and  maximum  velocity  difference,  which  is  the  basis  of  most  mixing  length 
theories.  This  implies  that  the  turbulent  properties  depend  on  the  axial 
distance,  are  constant  in  the  radial  direction,  and  are  In  local  equilibrium. 

One  might  expect  the  turbulence  to  lag  behind  the  equilibrium  value  where  the 
wake  is  growing  rapidly. 
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FIGURE  4  APPROXIMATE  TEMPERATURES  IN  WAKE  OF  MISSILE 


The  assumption  of  an  incompressible  core  becomes  less  valid  as  the 
point  at  which  the  core  and  inviscid  wakes  are  of  equal  width  is  approached. 

As  the  turbulent  core  reaches  the  edges  of  the  wake  the  velocity  difference 
in  the  turbulent  region  increases;  therefore,  the  turbulence  level  might  be 
expected  to  increase  in  this  region  of  increased  shear  and  cause  the  core  to 
grow  faster  than  predicted. 

The  wakes  computed  here,  while  only  approximate,  demonstrate  the 
turbulent  wake  is  much  shorter  than  the  laminar  wake.  The  physical  model 
and  the  value  of  con  be  evaluated  when  ballistic  range  data  are  available. 

Equation  (22)  predicts  a  constant  width  wake  at  largo  distances  from 

the  projectile  because  it  was  assumed  that  T^  «  T^.  Under  this  condition 

it  is  true  the  hot  turbulent  wake  will  grow  very  slowly;  however,  when 

T  T  the  wake  begins  to  grow  again  as  the  cube  root  of  axial  distance.  These 
00  ^ 

approximations,  except  the  assumption  of  local  equilibrium,  tend  to  produce 
wakes  which  grow  slower  and,  therefore,  are  longer  than  a  more  rigorous 
mixing  length  theory  would  predict.  Assuming  the  turbulence  to  be  in  local 
equilibrium  would  give  rise  to  shorter  wakes  than  a  lagging  turbulence  level. 
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NOMENCLATURE 


a,  b,  c,  X 

( 

D 

C 

P 

F 

h 

t 

L 


q 

r 

R 


T 


constants 
body  diameter 

specific  heat  at  constant  pressure 

drag 

enthalpy 

width  of  core 

width  of  Invlscld  wake 

heat  flux 

radial  distance 

body  radius 

U  I 

turbulent  Reynolds  number  =  — ^ 

€ 

temperature 


u  axial  velocity  In  core 

U  axial  velocity  In  Invlscld  flow  relative  to  atmosphere 

V  radial  velocity  In  Invlscld  flow 

^  axial  distance 

\  distance  at  which  i  -  L  behind  the  body 

^ ^  turbulent  viscosity 

^  turbulent  conductivity 

P  density 

T  shear 
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Subscripts 


a  value  at  r  =  o 

0  value  at  X  =  0,  V  =  o 

®  free  stream  value 
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